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We investigate theoretically and experimentally the statistical properties of dc current through an 
open quantum dot subject to ac excitation of a shape-defining gate. The symmetries of rectification 
current and photovoltaic current with respect to applied magnetic field are examined. Theory and 
experiment are found to be in good agreement throughout a broad range of frequency and ac power, 
ranging from adiabatic to nonadiabatic regimes. 



Transport in mesoscopic systems subject to time- 
varying fields combines elements of non-equilibrium 
physics and quantum chaos. This combination extends 
the scope of mesoscopic physics and is likely to be im- 
portant in quantum information processing, where fast 
gating and quantum coherence are both required. Of 
particular importance is the ability to control external 
fields applied to the mesoscopic system and to distin- 
guish effects of these fields on quantum dynamics of the 
system. For example, two distinct contributions to direct 
current through an open quantum dot due to an oscillat- 
ing perturbation have been identified P, 0] and observed 
experimentally in Ref. 

In this Letter, we investigate the statistical proper- 
ties of dc currents resulting from an applied ac electric 
field over a wide range of excitation frequencies, paying 
particular attention to the presence or absence of sym- 
metry with respect to magnetic field in various regimes. 
Theoretical analysis is based on recently developed time- 
dependent random matrix theory 0, 0| . Experiments use 
a gate-defined GaAs quantum dot subject to ac excita- 
tion of a gate at MHz to GHz frequencies. At low exci- 
tation frequencies lu <^t^^ (rd is the electron dwell time 
in the dot) the present theoretical results are consistent 
with those obtained by adiabatic approximations P, . 
However, the analysis is applicable over a wider range of 
frequencies uj ^ i?T, where i?T — ^/tctoss is the Thou- 
less energy and Tcross is the electron crossing time of the 
dot. At higher frequencies u ^ E^:, the system may be 
studied by methods developed for bulk conductors 0, Q . 

Three distinct contributions to dc current through the 
dot can be identified, resulting from: i) an applied dc 
bias; ii) an ac bias at the excitation frequencyfi.e., rec- 
tification effects 0); Hi) photovoltaic effects We 
restrict our attention to one-parameter excitation, not- 
ing that while in the adiabatic regime one- and two- 
parameter excitations affect the system differently, be- 
yond the adiabatic regime, w > r^"^, the differences dis- 
appear j^. 

The Hamiltonian of electrons in the dot in the pres- 
ence of a magnetic flux $ is represented by a Hermi- 
tian M X M matrix H{t) = 7i$ -I- Vcosut, with the 



time independent part being a random realization 
of a matrix from a Gaussian unitary ensemble with the 
mean level spacing 6i , and V being a matrix from a Gaus- 
sian orthogonal ensemble characterized by the strength 
Co = ttTtV^/M^Si and MSi Et The parame- 

ter Co determines the energy displacement of an electron 
state due to the applied perturbation V. The contact 
between the left (right) lead and the dot contains A'l 
(A^i) open channels, we enumerate channels, a, in the 
left (a = 1 . . . iVi) and the right {a = Ni + 1 . . . N^h) con- 
tacts, A^ch = A^i -I- Ni-. The corresponding experimental 
setup is shown in Fig. ^ 

The dc current /* through the dot is determined by 
the scattering matrix [Sg,{t,t')]ai3, see Ref. 0: 
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Here 5t,t' — S{t — t') and 
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is the distribution function of electrons in channel a at 
temperature T and voltage Va{t). At sufhciently low 
frequencies uj ^ Ec/h [Ec is the dot charging energy) 
Va{t) is simply related to the bias V{t) across the dot: 
Va{t) — KaotVit)- Elements of the diagonal matrix A are 
kaa = N^/Nch for 1 < a < Ai, and Aaa = -Ni/N^h for 
Ni<a< A^ch- 

We consider the bias V{t) across the dot in the form 
y{t) — Vb + Kj cos(cjt + (pi)- The dc current through the 
dot to first order in dc bias Vq and ac bias is 
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where the first term represents the photovoltaic current 
J*^ = P{Va = 0), see Eqs. ((TJ and 01. The second and 
third terms in Eq. Q represent the contributions to the 
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current due to dc bias Vq and ac bias Kj, respectively: localization (4>o = hc/e) 



Go-JG*(i) 
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where Gq ~ N\N,:/Nch is the classical conductance and 
SGk{t) — J^^^'^ SGk{t)ujdt/2TT stands for time averaging 
of the "instantaneous conductance" {k = 0, 1) 

5Gt{t) = J drFkir) J dO cos{k{oje + ^u)) (5) 



^°(-^- sinh(.JrA) '^^(-) 



27rfcBrsin(tJT/2) 
hio svcAiij^kT^TT jK) 



We observe that 5G\{t) = SGo{t) cos{ujt + ipi) in the 
adiabatic limit flu; <^ max{NchSi, ksT} considered in 
Refs. 
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FIG. 1: Symmetry factor Sk = ('JGfc SG~^)/({5Gff) as a func- 
tion of frequency w for k = (upper panel) and k = 1 (lower panel) 
at two values of temperature T and power Co of the ac excitation. 
Inset: Micrograph of device and schematic picture of applied volt- 
ages. 

Below we study the variance of the photovoltaic cur- 
rent and the conductances and gf with respect 

to random realizations of the Hamiltonian 7i$. Follow- 
ing Refs. 0, 01 , we find in the limit iVch ^ 1 and at 
magnetic fields $ ^ ^^yj Nch/M destroying the weak 
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and (/ph*/ph*) = 0, see Ref. In Eqs. © and the 
angle brackets (...) stand for the averaging with respect 
to realizations of 7i$. Functions B^^^.^ and Bf^,.^ de- 
scribe the distribution function 0| of electrons in the 
dot in the presence of time-dependent electric fields and 



kernels K 
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t,t':9- 



describe the evolution of elec- 

t,t';T 

and Bf^,.^ and the kernels K^^,.g^ contain the diffu- 

son I?(ti, t2,T) = cxp ^— J^^ r(T, t)dtj or the Coopcron 

Cin,T2,t) = exp(-i/;;r(T,t)dT). Here, r(r,t) = 

7c + 7^5 + 4Go sin^ sin^(cjT/2), where 7^ — SiNch/^n 
is the electron escape rate and is the electron phase 
relaxation rate due to inelastic processes. 

In the experiment, the ac bias Vu, results from capaci- 
tive coupling between the leads and the gate on which the 
ac voltage is applied (see the inset in Fig. 1). Therefore 
is proportional to the amplitude of the ac voltage at 
the gate. Assuming that Co is linear in the applied power 
to the gates, we write = a^y^Co/jc- The coefficient 
tto) has units of voltage and is independent of realizations 
of the quantum dot (we disregard fiuctuations of V over 
different realizations of 7i$). Therefore, the correlation 
functions of the rectification current /* = 5* are de- 
termined by the correlators of 6Gf : 

2 Go 



7o 



(5G+* (5G^ 



(8) 



We also notice that in the limit iVch ^ 1 the correlation 
function of the photovoltaic current and the rectifi- 
cation current if vanishes 0. 

First we use Eqs. JT)) and (|HJl to analyze the magnetic 
field symmetry of the rectification current /*. Although 
in the adiabatic limit uj 7c/7i the rectification current 
if is symmetric with respect to magnetic field inversion 
(<& ~* — *&), at higher frequencies lu ^ 7c/^ the symmetry 
of /* is suppressed. Indeed, the magnetic field symmetry 
is related to the time inversion symmetry. For a harmonic 
field at frequency lu, the time- inversion symmetry holds 
only on time scales much smaller than l/uj. Transport 
through the system is determined by times of the order 
of h/jc and consequently the magnetic field symmetry 
of the rectification current breaks if hui ^ 7c. We plot 

the ratio = (scf^ SG^)/((SGff) as a function of 
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FIG. 2: Upper panel: Variance of the conductance as a function of 
the ac excitation power P = PoCo/'fe at lu/27t = 5.56 GHz (hui = 
7.27o) and the theoretical result of Eq. Q with A: = 0. We use Pq = 
9 X 10~* W and k^T = 5. 470. Lower panel: Symmetric (O) S'lid 
antisymmetric (v) current correlators as a function of ac excitation 
strength Co- Solid line shows variance of the photovoltaic current 
Eq. 13 with parameters fixed by the fit in the upper panel. The 
dashed and dotted lines show the symmetric and antisymmetric 
correlators of the rectification current Eq. |2J with = O.A5Tiu)/e 
and fi = 0. 



fiLij 1 7o in Fig. ^ "5*1 = 1 represents the symmetric current 
/* oc (5Gf with respect to magnetic field inversion. In 
the adiabatic regime this symmetry originates from the 
Onsager symmetry |1SJ of the dc conductivity, see Eq. © 
and Refs. the frequency increases, vanishes, 

signalhng the suppression of the magnetic field symme- 
try. Therefore, the absence of magnetic field symmetry 
no longer serves as a distinct feature of the photovoltaic 
current which allows one to distinguish and the 
rectification current /*. 

We notice that the magnetic field symmetry of the dc 
conductance i^G* is more sturdy than the symmetry of 
the rectification current, see Fig.^ Particularly, at tem- 
peratures k-^T > 7e, dc conductance SGq is nearly sym- 
metric at frequencies Hlu ^ k^T, since the dc correla- 
tion function is determined by processes on a time scale 
Ti/k^T. The symmetry is not fully suppressed even at 
uj » k-QT/Ti] the suppression depends on Cq/^c- 

We apply Eqs. ® - © to the analysis of the experi- 
ment 3|. The quantum dot used in the experiment has 
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FIG. 3: Symmetric (Q) and antisymmetric (v) current correlators 
at uj/2-K = 2.4 GHz as a function of power P = PoCo/'fe with 
Po = 2.5 X 10~^ W. Solid line shows variance of the photovoltaic 
current Eq. JHJ at temperature T = 5.47o/fcB and frequency oj = 
3.l7c/?i. The dashed and dotted lines show the symmetric and 
antisymmetric correlators of the rectification current Eq. jSJ with 
a^j = 4.7hui/e and ipi = 0. 



an area A ^ 0.7 /im^. Relevant energy scales are the 
Thouless energy « 160 /leV and the mean level spac- 
ing 61 = TTjc ~ 10 /xeV. The measurements were per- 
formed at the base electron temperature T w 200 mK 
(/cbT = 5.470). From the size of conductance fiuctua- 
tions without ac fiuctuation on the gate, we estimate the 
dephasing rate 7,^ « 0.27c (see fl^ for details) and thus 
disregard it in our quantitative analysis. 

In the upper panel of Fig. |21 we show the variance of 
the conductance as a function of the incident power P 
at w/27r = 5.56 GHz {huj/^e « 7.2). We also plot the 
variance of the conductance calculated from Eq. ^ {k — 
0) at temperature T = bAje/kB- Assuming that the ratio 
Co /7c is proportional to the power P of the ac excitation 
applied to the gate, i.e. Cq/jc = P/Pq, we rescale P 
to obtain the best fit of the experimental points by the 
curve of Eq. (jT)). We find Pq = 9 x 10"*^ W. 

In the lower panel of Fig. |21 we show the correla- 
tors (J"'"*/^*) of the measured current. Although the 
traces of the magnetic field sweeps look quite asymmetric 
for the measured current, the antisymmetric correlator 
(/"'"*/"*) is not significantly smaller than its symmetric 
counterpart. We notice however, that if the averaging 
is performed over n realizations, the measured correlator 
(/+*/-*) can be estimated as (/+*/+*)/ (n - 50 in 
the experiment). 

We plot the variance of the photovoltaic current, 
Eq. ©, as a function of Co/^e for /cbP — 5.47e and 
huj ~ 7.27c (to facilitate numerics, wc used the approx- 
imation hto 3> 7c). We emphasize that the horizontal 
shift between the data points and curve is fixed by the 
fit in the upper panel for the dc conductance and there 
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is no fitting parameters for the variance of the current 
(along vertical axis). At Cq ^ 7e, the variance of the 
measured current changes quadratically in Co/70, consis- 
tent with quadratic dependence on Cq of the theoretical 
curves for (/p^) and therefore our assumption that Co is 
proportional to the power of the ac excitation is justified. 
At Co > 7o the variance of the measured current starts 
saturating. This saturation is expected for large power 
asymptote of the photovoltaic current due to spreading 
of the distribution function of electrons in the dot pj- 
Some deviation between the experimental points and the 
theoretical curve is expected due to the approximation 
A^ch ^ 1 used for derivation of Eq. © {N^h = 2 in the 
experiment). 

For illustration, we also plot the correlation functions 
of the rectification current, using Eq. JSJ for ipi ^ and 
au: = OAbhiv/e. For the rectification current, the satura- 
tion at large power is not expected: according to Fig. |2 
(J+*/±*) cx (Co/7o)° with a « 0.6. 

We similarly discuss the data for uj/2Tr = 2.4 GHz. 
Performing the fit of the experimental values of the con- 
ductance fluctuations and the result of Eq. ([TJl with k = 
and temperature T = BAje/kB, we find the relation be- 
tween the strength of the perturbation Co and the power 
P = PoCo/7c with Po = 2.5 X 10-^ W. 

In Fig. |31 we show the symmetric and antisymmetric 
current correlators for 2.4 GHz. For comparison we plot 
by a solid line the variance of the photovoltaic current 
/pij, calculated from Eq. © at a; = 3.l7c. We observe 
that the fluctuations of the measured current significantly 
exceed (by a factor ~ 100) the expected magnitude for 
the photovoltaic current, and therefore are likely due to 
the rectification of the bias across the dot. The low power 
data can be fitted by Eq. ^ with a^^ = A.Thuj/e. 

The above choice for a^^ > {S-w, k^T} je limits the ap- 
plicability of the linear expansion Eq. ||2Jl to small powers 
of the ac excitation Co, such that Co/7c ^ [Uuj / ea^Y . 
The higher order corrections in the bias VL, do not restore 
magnetic field symmetry, which is in apparent contradic- 
tion to the observed symmetry of the measured current at 
larger powers (at Co/7e ^ 1 in Fig. We attribute the 
restoration of magnetic field symmetry to dephasing due 
to dot heating by the dissipative current. Increasing the 
power P at fixed oj drives the system into the adiabatic 
regime since the heating makes the ratio ?iw/(7c -t- 7ip) 
decrease. As shown already in Fig. ^ the rectification 
current is symmetric in the adiabatic regime. The as- 
sumption that 7;p increases as power P increases is con- 
sistent with the observed change of the correlation field 
for the current fluctuations, see Fig. 4 in Ref. 3]. 

In summary, we studied ensemble fluctuations of dc 
current through an open quantum dot subject to oscil- 
lating perturbation. We showed that as frequency of the 
perturbation increases, magnetic fleld symmetry of the 
current disappears, regardless of the mechanism of the 



current generation. We demonstrated that the power be- 
havior of the current fluctuations is an important tool to 
distinguish effects of an ac excitation on dc current. 
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